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ABSTRACT 

We perform cosmological hydrodynamic simulations to determine to what extent galaxies 
lose their gas due to photoheating from an ionizing background. We find that the characteristic 
mass at which haloes on average have lost half of their baryons is M c ~ 6.5 x 10 9 hr 1 M at 
2 = 0, which corresponds to a circular velocity of 25 km s _1 . This is significantly lower than 
the filtering mass obtained by the linear theory, which is often used in semianalytical mod- 
els of galaxy formation. We demonstrate it is the gas temperature at the virial radius which 
determines whether a halo can accrete gas. A simple model that follows the merger history 
of the dark matter progenitors, and where gas accretion is not allowed when this temperature 
is higher than the virial temperature of the halo, reproduces the results from the simulation 
remarkably well. This model can be applied to any reionization history, and is easy to incor- 
porate in semianalytical models. 
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1 INTRODUCTION 

A UV-background produced by quasars and stars acts to quench 
star formation in small galaxies by photoheating their gas, which 
gets too hot to be confined in th eir shallow potentia l wells . 
This was originally pointed out by iDoroshkevich et al] d 1967b . 
and investi gated in the context of t he cold dark matter (CDM) 
model by ICouchman & Reesl {T986). Depending on the spectral 
shape of the ionizing background, the typical virial temperature 
below which haloes lose their gas is around 10 4 K. Photoheating 
associated with reionization therefore inhibits star formation in 
dwarf galaxi es and affects the faint end of the galaxy luminosity 
function (e.g. lEfstathiou]|l992l) . 



ISusa & Umemural d2004allbh . 

However there is still debate about the value of the char- 
acteristic mass, M c , below which galaxies are strongly affected 
by photoionization. iGnedinl <|2000) argued that Af c = Mf, the 
filtering mass that corresponds to the scale over which baryonic 
perturbations are smooth e d in linear perturbation theory (see 
Appendix FBI. Hoeftetal. (2006) used simulations to argue that 
M c -C Mp, in particular at low redshift. They argued that 
M c follows from considering the equilibrium temperature, T oq , 
between photoheating and radiative cooling at a characteristic 
overdensity of A ~ 1000; M c haloes have a virial temperature 



T cq (A : 



1000). 



Since iReesI dl986h argued that in dark matter haloes with 
circular velocities around 30 km s _1 gas can be confined in a 
stable fashion, with radiative cooling balancing photoheating, 
many attempts to quantify the effects of an ionizing background 
on galaxy fo rmation have been made using semianalytical calcu- 



lations (e.g. Babul & Rees 1992; Efstathiou 1992; Shapiro et al 
19941; iNagashima et al] 1 19991 : iBullock et all l2000l : IBenson et al 
2002al ltt" Somerville 2002J), spherically symmetr ic 
dThoul & Weinberg! 11996k iKitavama etatl |2000|) . and three 



spherically symmetr ic simulations 
atavama et all l2000t) . 
nydrodynamic simulation ; 
19961: iNavarro & Steinmetzf [l997l; IWeinberg et ail [ 



dimension al cosmological hydrodynamic simulation s (Ouinn et al 

" 1991 iGnedin 



2000). The effects of self-shielding have been investigated us- 
ing three-dimensional radiative hydrodynamic simulations by 
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The relation M c = Mf proposed by IGnedinl is often 
used in semianalytic modelling for dwarf galaxy formation in 
order to describe the quenc hi ng of star formatio n dBullock et al.l 
l2000l : IBenson et al] l2002j|bl: ISomervillel |2002|) . An important 
application of this relation is to estimate whether reionization 
can explain the apparent dearth of satellite galaxies in the Milky 
Way, a s compared with th e high abundan ce of dark matter sub- 
haloes dMooreetai]|l999l) . Interestingly, ISomervillel d2002l) and 
Nagashima & Okamot ol d2006t) argued that models underestimate 
the number of dwarf satellite galaxies in the Local Group if the 
effect is as strong as inferred from the filtering mass. 

Given the uncertainty in M c , we were motivated to perform 
high-resolution hydrodynamic simulations of galaxy formation in 
a ACDM universe with a time-evolving UV-background, and to 
measure the evolution of the baryon fraction as function of halo 
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mass. We compare our results to those of lGnedir] and lHoeft et all 
and formulate an intuitive and simple model that reproduces our 
simulations very well. The model can be incorporated easily into 
any semianalytic model of galaxy formation. 

This paper is organised as follows. Details of our simulations 
are described in Section 2. Section 3 contains our numerical re- 
sults which focus on the baryon fraction of the simulated haloes, 
and the characteristic mass at which the baryon fraction is on aver- 
age half the cosmic mean. We present a simple physical model that 
describes the evolution of the baryon fraction of individual haloes 
and compare it to the simulations in Section 4. Discussion and a 
summary follow in Sections 5 and 6. 



2.2 Runs 

We use a modified version of the PM-TreeSPH code GADGET2 



2 SIMULATIONS 

2.1 Background cosmology and radiative processes 

The physics of how low-mass haloes fail to accrete gas, or 
lose it after the gas is being ionized, is largely independent 
of cosmology, but for definiteness we assume a geometrically 
flat, low density CDM universe (ACDM) with cosmological 
parameters: fi = 0.25, «a = A. /(3Hg) = 0.75, ft b = 0.045, 
h = #o/(100 km s" 1 Mpc" 1 ) = 0.73, and <j 8 = 0.9. Symbols 
have their usual meaning. 

The radiative processes that we take into account are Compton 
cooling off the microwave background, thermal Bremsstrahlung 
cooling, and line cooling and photoionization heating from Hydro- 
gen and Helium i n the presence of a n imposed ionizing background 
as computed by lHaardt & Madaul d200 lh ■ Although our cooling 
routine handles cooling and photoheating for many species, we 
assume primordial gas throughout this paper for simplicity. The 
routine was developed for a different project and will be described 
elsewherfj Bri efly it employs tab ulated reaction rates evaluated 
with CLOUDY dFerland et alJI 199*3) . assuming ionization equilib- 
rium. We assume Hydrogen and Helium I reionization occurs at 
z = 9, and Hell reionization at z = 3.5. 



We increase the thermal energy by 2 eV per atom both 
during H and Hell reioni zation to mimic r adiativ e transfer and 
non-equilibrium effects jAbel & Haehneitl 1 19991) . This brings 
the evolution of To and 7eos that characterise the low-density 
temperature-density rel ation, T(p) = Tp ( p/(p)) 7EOS ~ 1 closer 
to the measurements of ISchave et alj dl99S ) n. The extra energy 
is distributed over redshifts with a half-Gaussian distribution 
with a dispersion a z = 0.0001 from z — 9 for Hydrogen 
reionization and with a Gaussian distribution centred on z = 3.5 
and with a z — 0.5 for Hell reionization. While our numerical 
results depend on the assumed ionizing background, our model, de- 
scribed in Section[4] can be applied to any assumed thermal history. 



Springell 12005 ), the successor of the TreeSPH code GADGET 



Springel et alj 



200 lh . Hydrodynamics is treated with Smoothed 



Partic le Hydrodynamics (SPH; lLucvlll977l ; lGingold & Monaghan 
1 1977b . employing a 'conservative en tropy' formulation that man i- 
festly conserves energy and entropy (Springel & Hernquist 2002). 

We have performed a series of runs which differ in numerical 
resolution and box size (TableQJ. We also use a run with a different 
imposed ionizing background. The 'reference' simulation uses 
a periodic cubic box of comoving size L = 4 h~ x Mpc, and 
equal numbers of TV" = 256 3 dark matter and SPH particles. The 
'high-resolution' and 'low-resolution' simulations are used to 
gauge numerical convergence. The 'high-resolution' simulation 
used zoomed initial conditions, where a region of size Lh = 2h~ 1 
Mpc inside the 'reference' box is simulated with eight times better 
mass resolution. The region outside the high-resolution region is 
populated with dark matter particles only, at coarse resolution. 
Linear waves common to both simulations have identical phase 
and amplitude. The 'low-resolution' simulation uses a 8 h,- 1 Mpc 
comoving cube, with uniform mass resolution, eight times coarser 
than the 'reference' simulation. Simulation 'reference-lbox' has 
identical resolution to 'reference', using zoomed initial conditions 
in a larger box. 

We convert SPH particles into collisionless ('star') par- 
ticles in high density regions (those with gas overdensity 
A ee p/(p) > 1000 and physical Hydrogen number density 
7iH > 5 cm -3 ) in a time given by the local dynamical time, 
tdyn = (47rGp) _1,/2 . This decreases the computational cost of the 
simulations significantl}[fl and we have verified it does not change 
the results as compared to simulations without star formation. 

The evolution of gas in the temperature-density plane is 
shown in Fig.[T] Before reionization, low density gas lies along the 
adiabat T oc p 7-1 , with 7 = 5/3, except where shocked, whereas 
cooling introduces a sharp maximum temperature of T ~ 10 4 K at 
higher density. Bigger simulation boxes will contain more massive 
haloes which do contain shocked gas above 10 4 K. Heat input 
during reionization makes the gas nearly isothermal at T ~ 10 4 K. 
At later times gas is above a minimum temperature T m i n (p, z), set 
by the balance between photoheating and adiabatic cooling at low 
density, and photoheating and radiative c ooling at higher density 
dHui& Gnedir]|l99l iTheuns et alj|l998h . Shocked gas forms the 
plume above T m j n (p, z = 5) seen in Fig.Q] 

The equilibrium temperature T eq (p, z) where photoheating 
balances radiative cooling is indicated by crosses for overdensities 
A = 200 and 1000, whereas we also include adiabatic cooling for 
the cross at A = 1. As expected these points fall on the T , m i n (p, z) 
curves both for z = 8.71 (just after reionization) and the later 
redshift z = 5. 



1 We would like to thank our colleagues J. Schaye, C. Dalla Vecchia and 
R. Wiersma for allowing us to use these rates. 

2 Throughout this paper, p and (p) denote physical gas density and the 
mean baryon density in the universe, respectively, unless otherwise stated. 



3 We do not include stellar feedback, since we concentrate here on radiative 
feedback only. 
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Table 1. Numerical parameters of the performed runs. Simulations use periodic boxes with the indicated box size, 'zoomed' simulations evolve part of the 
simulation box at higher resolution. mgpH and mju denote the particle masses of SPH and (high-resolution) dark matter particles, respectively, e is the 
comoving gravitational softening length in terms of the equivalent Plummer softening. The gravitational force obeys the exact r~ 2 law at r > 2.8e. 



Simulation 
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Zoom 


-^SPH 


"SPH 
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(IT 1 Mpc) 






(10 4 /i- 1 M Q ) 


(10 4 fc- 1 M Q ) 


(h' 1 kpc) 


reference 


4 


No 


256 3 


4.76 


21.7 


0.5 


high-resolution 


4 


Yes 


256 3 


0.596 


2.71 


0.25 


low-resolution 


8 


No 


256 3 


38.1 


174 


1 


reference-lbox 


8 


Yes 


256 3 


4.76 


21.7 


0.5 




Figure 1. Distribution of the particles in the density-temperature phase diagram for the 'reference' simulation, colour coded according to the fraction 
of particles at given p and T, at redshifts z = 9.31 (before reionization), z = 8.71 (just after reionization) and z = 5, as indicated in the panels. 
Adiabatic cooling, impulsive heating during reionization, photoheating, radiative cooling and shocks from structure formation determine the distribution and 
its evolution. Plus signs indicate temperatures obtained assuming gas evolves at constant overdensity. This reproduces T(z) well for the indicated overdensities 
of A = 1, 200, and 1000. 



3 ANALYSIS OF SIMULATIONS 



3.1 Halo identification 



haloes with at least 100 dark matter particles, and discard haloes 
contaminated by boundary particles in the zoomed simulations. 



We identify haloes and the baryons associated with them as follows. 
First we identify overdense region s using the friends-of-friends 
(FoF) algorithm dDavisetal.ll 19851) on t he dark matter particles 
only, choosing the linking according to [okamoto & Nag ashim3 
J2003h . This picks-out overdensities according to the spherical top- 
hat model appropriate for the assumed cosmological parameters 
dEke et al.lll996h . 

We refine the location of the centre of each given halo, by 
determining the centre of mass of all particles w ithin a sphere 
of giv en radius R, and progressively shrinking R jLacev & Cold 
1 19941) . The location of the density peak is identified as the centre 
of the halo. Finally we determine the virial radius of the halo, 
R v ii, such that the mean dark matter density within R v [ r equals 
the virial density of dark matter component. 

Finally, SPH and star particles within R v i r from the halo 
centre are assigned to that halo. We restrict the analysis below to 



3.2 Baryon fraction of simulated haloes 

Before reionization the baryon fraction /b = Mb /M of simulated 
haloes scatters around the cosmic mean, (/b) = fib/^o> as ex- 
pected (Fig[2] the drop toward lower masses is a resolution effect). 
After reionization /b drops sharply below the characteristic mass 
M c (z)', haloes with M ^ M c (z) have too shallow potential wells 
to hold on to their photoheated gas. The drop of ft, from the cosmic 
mean to just of few per cent occurs over approximately a decade 
in total mass. The characteristic mass M c where /b = (A>)/2 
is the same for 'reference' and 'high-resolution' simulations 
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Figure 2. Fraction of baryonic mass as function of total mass, M, for simulated haloes at three redshifts z for all haloes (top row) and 'isolated haloes' (bottom 
row). The cosmological baryon fraction Qf,/ Qq m 0. 18 is indicated as the dotted horizontal line. Plus signs and diamonds indicate haloes from the 'reference' 
and 'high-resolution' simulations, respectively. Their comparison shows that the drop in M b /M below M ~ 10 8 Mq before reionization (z = 9 in 
these simulations) is an artifact of numerical resolution; at higher masses or after reionization both simulations give similar values, demonstrating numerical 
convergence. The baryon fraction M^/M falls from the cosmic mean to just a few per cent over about a decade in mass below a critical value M c which 
depends on z. Halo to halo scatter is much reduced when considering only isolated haloes (compare top with bottom panels), which removes tidally stripped 
haloes. Full lines indicate fits using equation (TJ to the 'reference' model; these fit the simulated haloes very well. The dashed lines in the bottom panel repeat 
the fits to the top row; the critical mass in isolated haloes is slightly higher than for all haloes. 



demonstrating numerical convergence- 
Some haloes have fb S> {fb)- By tracking their merger 
histories we find that most of these were subhaloes (i.e. part of 
a larger halo) at earlier times, which have been tidally s t ripped 
from their hosts leaving baryon-rich cores. iLudlowet all (2008) 
found that many apparently isolated low-mass haloes are subhaloes 
ejected by multiple-body interactions. This phenomenon of tidal 
stripping may have occurred for haloes with /b ^ (fb) as well; 
therefore the baryon fraction of such low-mass haloes may be 
artificially high. 

To reduce the contamination from tidally stripped haloes we 
define isolated haloes as those that lie outside 6 R v i r of more 
massive haloes. The baryon fraction of isolated haloes shows much 
less scatter as function of halo mass (Fig. [2] bottom panels), and 



4 We stop the 'high-resolution' simulation at z = 5 because of its small 
volume: the number of haloes with M > M c becomes too small to de- 
termine M c . We stop the 'reference' simulation at z = 2.09 for the same 
reason. 



we restrict our analysis to isolated haloe^f]. 

We use the following function proposed by GnedijJ ( feOOOh to 
describe how the baryon fraction depends on mass and redshift: 

/ b (M,,)H/ b ){l + (2^-l)(^)" a } ' (1) 

in terms of the cosmic mean baryon fraction, (fb)- This function 
has two fitting parameters: a and M c (z). The baryon fraction 
of haloes with M 2> Af c equals the cosmic mean while that of 
haloes with M <C M c goes to zero oc (M/M c ) 3 . The parameter 
a controls how rapidly fb drops for low mass haloes; a value of 
a = 2 fits the si mulations well (Fi g. [2j, and is also consistent 
with the results of iHoeft et al.l (j2006). Haloes with mass equal to 
the characteristic mass, M = M c have lost half their baryons. As 
expected, the value of M c for isolated haloes is slightly higher 
than for the full sample of haloes (by 27 per cent at z — 5 and 50 
per cent at z = 2.09 for the 'reference' simulation; see Fig. |2j. 



5 This decreases the number of haloes from 3618 to 1185 for the 'refer- 
ence' simulation at z = 2.09. 
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Figure 3. Top panel: Evolution of the characteristic mass M c with red- 
shift z. Different symbols refer to simulations with the same assumed 
ionizing background but different resolution and/or size of the simulation 
box( 'reference': circles; 'high-resolution': stars; 'low-resolution': squares; 
'reference-lbox': triangles). When haloes of M = M c contain fewer than 
10 3 dark matter particles the simulation is shown with open symbols. Over- 
lap of the different filled symbols demonstrates numerical convergence. M c 
rises from ~ 10 7 h~ 1 M Q just after reionization to ~ 6.49 X 10 9 h _1 M Q 
at z = 0. Lines indicate the virial masses corresponding to virial temper- 
atures of T vir = T cq (A = l/3A vir ) (dotted) and T cq (A = 1000) 
(dashed); the redshift dependence of these lines is quite different from that 
of M c , demonstrating a single temperature cut is not sufficient to com- 
pute the characteristic mass. It is necessary to follow the merger history of 
the haloes. Lower panel: corresponding evolution of the virial temperature 
(left hand scale) and circular velocity (right hand scale). M c has been con- 
verted to T v i r using Eq. (4j. The characteristic circular velocity V C (M C ) 
rises rapidly to f» 10 km s~ x following reionization, then continues to 
increase more gradually reaching ?s 20 km s —1 at z = 2, and remains 
roughly constant thereafter. The equilibrium temperatures at two overden- 
sities (A = 1/3 A v ; r : dotted; A = 1000: dashed) are also shown; they 
are nearly constant, whereas the virial temperature of haloes with mass M c 
increases rapidly especially at early times. 



The characteristic mass increases from M c [z] 

t9 



10 7 h^M® 

just after reionization to M c as 6.49 x 10 y h~ 1 M Q at z = 
(Fig. [3](- The 'high-resolution' simulation has a smaller value of 
M c by 20 per cent as compared to the 'reference' simulation, 
at z = 5. This may simply reflect the small number of massive 
haloes in the 'high-resolution' simulation (see the lower middle 
panel of Fig. [2}. Haloes with fewer than 10 3 dark matter particles 



may suffer from two-body heating due to massive dark particles 
( Stein metz & White] 1 19971) ; they are shown as open symbols. 
Values of M c for these less well-resolved haloes are only slightly 
higher than of haloes resolved with more particles. Overall though, 
simulations that differ in numerical resolution and/or box size give 
similar results, demonstrating numerical convergence. This gives 
confidence in the value of M c (z = 0) » 6.48 x 10 9 h~ 1 M Q 
from the 'low-resolution' simulation. The a mazing agreement w ith 
the value of 6.5 x 10 9 h~ 1 M Q obtained bv lHoeft et all J2006h for 
simulations of voids may be slightly fortuitous. 

We convert halo mass M to virial temperature T v i r using 



1 /im p 2 

M3 ) 



2 k B 



(2) 



where 



(3) 



/ GM\ 5 

v - = («-) 

is the circular velocity of the halo at the virial radius. While we 
compute the mean molecular weight p self-consistently in our sim- 
ulations, here we use p = 0.59 to compute T v i r for simplicity. In 
terms of the overdensity A v i r within the virial radius R v i r , Tvir 
depends on halo mass and redshift as 



T vir (M,z) 



1 pm p 
2~feT 



A v i r fin 



(l + z)(GMH ) 



(4) 



The value of T v i r at the critical mass M c (z) rises rapidly to ~ 
5000 K following reionization, keeps increasing to ~ 10 4 K at 
z = 2, then increases much more slowly to ~ 2.2 x 10 4 K (V c ~ 
25 km s~ 1 )atz = (Fig. [3] bottom panel). Our z — value 
is much lower than the temperature below which galaxy formation 
is assumed to be strongly suppres sed in semianalytica l models of 
galaxy formation (e.g. ~ 10 5 K in lBenson et alj|2002bl) . 



4 MODELLING THE SUPPRESSION OF THE BARYON 
FRACTION BY THE PHOTOIONIZING 
BACKGROUND 

A simple model of gas accretion compares the halo's virial 
temperature T v i r with the gas temperature T: if T > T v i r , there 
is no accretion. However the temperature T of the IGM depends 
on its density p (Fig.Q}: at which density should we evaluate T? 
Gnedin's filtering mass sug gests to use a densi ty close to the mean 
(see Appendix [B} whereas Hoe ft et ail d2006h used p = 10 3 (p) 
appropriate for significantly overdense gas inside haloes. However, 
if one considers the density profile of a halo, gas at the edge of 
the halo (at J? v ir) has a local density of order one-third of the 
virial gas density, p v - u /3. This suggests that gas will be prevented 
from accreting when T cq (p v ir/3) > T v i r . This temperature, 
Teq(pvir/3), is uniquely defined as a function of redshifjf] from 
cooling and heating functions used in the simulations. 

Simply identifying M c from the condition T Gq (p v i r /3) = T v i r 
does not work well, overestimating M c (z — 0) by a factor ~ 2 as 
compared to the simulations (dotted line in Fig. [3}. In addition, the 
redshift dependence of M c (z) is quite different as well. Clearly 
this model is too simple, as it ignores the baryon fraction of the 



Note that the equilibrium temperature also depends on gas metallicity. 
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Figure 4. Baryon fraction of isolated haloes calculated by our models compared to values from the SPH simulation for isolated haloes with more than 100 
dark matter particles and with / b > /™ m picked from the 'reference' simulation at z = 2.09. The simple model A (left hand panel) typically reproduces the 
baryon fraction obtained from the simulation within 50 per cent, whereas model B (right hand panel) improves the agreement especially for haloes having low 
baryon fraction, / b < (/ b )/2. The dashed lines indicate |/ b: modcl - / bj sim |// b , sim = 0.5. 



progenitors of a given halo, as well as the photoevaporation of 
gas from haloes during reionization. F or the same reaso ns, the 
agreement at z = between a model bv lHoeftet ail (2006) which 
compared the equilibrium temperature at high density (A = 1000) 
with T v i r (dashed line in Fig. [3} and the simulations is probably 
a coincidence. Next we describe an improved model that takes 
the merger history of haloes into account, as well as modelling 
photoevaporation and the decrease in accretion rate if gas it too hot. 

Before reionization we assume each halo has the cosmic 
baryon fraction, (/ b ). We begin to construct merger trees of FoF 
haloes with more than 32 dark matter particles once the Universe 
is reionized. We model photoevaporation during reionization 
by computing the equilibrium temperature at large overden- 
sity and comparing it to the halo's virial temperature: when 
r cq (A CV p) > Tvir the baryon fraction decreases oc exp(— t/t cvp ), 
where the evaporation time-scale t cvp = i? v j r /c s (A cvp ). The 
time t is counted either starting from reionization, or from the 
previous merger (whichever is later), R v i r is the virial radius 
and the sound speed c s is evaluated at T e q(A evp ). Since T cq (p) 
is nearly independent of p for high densities, the exact value of 
the overdensity is not important as long as A cvp ^> A v i r ; for 
reference we use A evp = 10 6 . 

When two haloes merge, we add the baryon masses of the 
progenitors, and we also allow the halo to 'accrete' baryons, but 
only when the temperature of the accreting gas is lower than the 
virial temperature: T eq (p v ir/3) ^ T V - 1T . Haloes that form after 
reionization at the resolution limit of the FoF are assumed to 
have the minimum baryon fraction^ /b™. This completes our 
model 'A'; it reproduces / b (M) as measured from the simulations 

7 We assume the baryon density in a halo is at least the mean baryon den- 



reasonable well (Fig[4] left hand panel), with a small tendency to 
under-predict / b for haloes with / b ~ 0.4 in the simulations. The 
Appendix describes a more involved model 'B' that does better for 
such haloes (FigO right hand panel). 

In detail, model A works as follows. Before reionization each 
halo has the maximum baryon mass given by the cosmic mean 

</b>: 

M b = M b max = A/dm . (5) 

1 - [Jb) 

The total mass of the halo is simply the sum of its baryonic and dark 
matter masses: M = M b + Mdm. After reionization the baryon 
mass of a halo that forms from the merger of several progenitors is 

M b = M b + M acc 

pr ° s / ft \ 
M b = Vexp --MA4, (6) 

\ £ev P / 

the sum of baryon masses of its progenitors, A/ b , and the accreted 
baryons, Af acc . The decrease in baryon fraction exp(— St/t evp ) is 
due to photoevaporation, with St the elapsed time since the previous 
merger. For a progenitor with T cq (A cvp ) < T v ir, we set St/t cvp = 
since evaporation does not take place in such a halo. The mass 
available for accretion is given by 

M acc = max(M b max - M b) 0) . (7) 

When the temperature of the accreting gas, T acc = T , eq (p v ir/3) is 
greater than the virial temperature of the halo, T v i r , the halo does 
not accrete any gas: A/ acc = . 



sity in the universe. Thus / b nin = (p)/(p°^ + (p)) < W~ 2 . This gives 
a good estimate of the baryon fraction of low-mass haloes with M M c 
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Figure 5. Evolution of the characteristic mass, M c , below which haloes 
half lost half of their baryons, comparing numerical simulations (symbols) 
with the simple models A and B described in the text. These models use 
the dark matter merger trees with simple recipes for when haloes lose and 
accrete gas. Both models do well, with M c (z) from model B nearly in- 
distinguishable from the full simulation. The middle panel also illustrates 
model B but neglecting photoevaporation during reionization: this leads to 
a significant underestimate in M c . 
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Figure 6. Upper panel: Temperature evolution of gas a t several overden- 
sities indicated in the legend, for simulations with the Haardt & Madau 
1 2001 ) photoheating rate (thin lines), and an experiment with increased pho- 
toheating at high density (thick lines). Lower panel: Filtering mass, Mp, 
computed from (T)v and characteristic mass, M c , for the original photo- 
heating rate (dotted line and circles) and modified photoheating rate (full 
line and diamonds). (Open circles indicate that haloes of total mass M c 
do not contain a sufficient number of dark matter particles for reliable mea- 
surements.) Modifying the photoheating rate does not change Mp, but does 
increase M c significantly. Model A with the modified heating rate (dashed 
line) reproduces this increase well. (Model B is not shown; it gives nearly 
identical results to model A). 



We can now predict the baryon fraction for all haloes selected 
from a simulation at a given redshift, and fit Eq. |QJ to the dis- 
tribution of fb(M, z); this yields the characteristic mass M c (z). 
This function is compared to that obtained from the simulations 
directly in Fig.[5j model A gets the simulated value to better than 
20 per cent whereas model B is almost indistinguishable from the 
simulation. We also show model B without taking photoevapora- 
tion into account; this leads to an underestimate of M c by more 
than a factor of two around z w 6, illustrating the importance of 
including this effect. 



5 DISCUSSION 

Our numer i cal re sults are in good agreement with those of 
iHoeft et alj J2006h who used the same simulation code (GAD- 
GET 2) but a different implem entation of photoionization and cool- 
ing faaardt & Madaul j 19961) UV-background with reionization at 
z = 6). The characteristic mass, M c , at which haloes lose half of 
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their baryons is much less than the value found bv lGnedinl J200oh : 
we find that reionization is not effective at suppressing gas accre- 
tion in haloes of the filtering mass Mf (Appendix |Bj, only haloes 
with M sC M c <C Mf are strongly affected. Our simulations sug- 
gest M c is mainly determined by the temperature of the accreting 
gas at a density of on e-third of the virial density, A ~ l/3A v i r , 
whilst [Gnedinl (120001 argued it is the temperature of the IGM at 
the mean densit)QTb ■ The values of To in his and our sim ulations 
are in fact quite similar. A much stronger effect found bv lGnedinl 
feOOd) was probably due to insufficient resolution (N. Gnedin, pri- 
vate communication). 

It is however important to note that radiative transfer may lead 
to spectral hardening and may preferentially heat gas in overdense 
regions. To mimic this effect, we ran a simulation where we artifi- 
cially increased the photoheating rate by a factor of 10, for gas at 
overdensity A > 10: 



Turned (A, z) 



H(A,z) 
WH(A,z) 



A < 10 
A > 10. 



(8) 



Ti(A, z) is the original rate at overdensity A and redshift z, and 
7i mo d(A, z) is the modified rate. This arbitrary modification heats 
gas at A = l/3A v j r to much higher temperatures yet does not 
affect To (Fig. [6]l. Since Tq is not affected, neither is the filtering 
mass. However since the accreting gas is hotter, our model predicts 
that M c will be increased. 



resolved the haloes whose mass is below the filtering mass, these 
simulations were stopped just after the reionization where the dis- 
crepancy between the filtering mass and the characteristic mass is 
small (see Appendix [B]l. We expect their model can be greatly im- 
proved by using our model of gas accretion and evaporation. 



6 CONCLUSIONS 

We have performed high-resolution cosmological hydrodynamical 
simulations with an imposed, spatially uniform but time dependent 
UV-background to investigate to what extent haloes lose baryons 
due to photoevaporation. Our numerica l results ar e con sistent 
with the simulations of void regions by iHoeft et alj (2006), but 
inconsistent with those found earlier by Gnedin (2000). The effect 
of a photoionizing background on galaxy formation is much 
weaker than previously assumed in s e mianalytic model s (e.g. 
iBullock et al.l200d : lBenson et al.l2002allbl:ISornervillei2002l) based 
on the filtering mass approach proposed by iGnedinl (2000) (see 
Appendix iBt, 



The volume-averaged temperature (T)v in the 'reference' 
simulation with the modified photoheating is shown in the upper 
panel of Fig. [6] and is as expected almost identical to that in the 
original 'reference' simulation. Consequently, the filtering mass 
obtained with the modified photoheating is indistinguishable from 
that in the 'reference' simulation even if we compute Mf from 
(T)v instead of To (lower panel). The characteristic mass is 
however significantly higher, as predicted by our model. 

This numerical experiment provides three important insights: 

• The mean IGM temperature is not important for the evolution 
of the baryon fraction of the haloes. Hence, any agreement between 
the characteristic and filtering mass must be a coincidence. 

• The characteristic mass is sensitive to the temperature at the 
overdensity around A = l/3A v j r , and therefore spectral harden- 
ing from radiative transfer effects can, in principle, increase the 
characteristic mass significantly. 

• Our model works well for any reionization histories as long as 
the gas temperature as a function of density and redshift is known. 

The modified heating rate used above probably significantly 
over estimates the effect of spectral hardening. Although radiative 
transfer effects may be important during reionization, their effect 
will lessen after the gas is ionized. Therefore we believe our results 
with the original heating rate are robust, at least for lower redshifts 
(z < 6). 

Interestingly, Ras era & Tevssiej {2006) found a good agree- 
ment between their simulations and their model on the baryon bud- 
get, which was based on the filtering mass. The reason for this 
agreement is probably that the minimum halo mass above which 
baryons are able to fall into the dark matter potential well was con- 
trolled by the numerical resolution rather than the filtering mass as 
they discussed. Although their highest resolution simulations well 



Gnedin used the volume-averaged temperature (T) v instead of To 



We have presented a novel model which uses the merger 
tree of each halo, and very simple physical arguments about gas 
retention and accretion, to predict how much the baryonic mass is 
in haloes of given total mass. The model is in excellent agreement 
with the simulations. Gas accretion is modelled very accurately 
if one simply assumes that gas at the outskirts of the halo, at an 
overdensity of a third of the virial gas density, cannot accrete if it is 
hotter than the halo's virial temperature. A slightly more involved 
model does even better in predicting baryon fractions as function 
of halo mass and merger history, as compared to simulations. 
These models allow us to predict the baryon fraction of individual 
haloes and hence the characteristic mass M c (z), at which haloes 
lose half of their baryons, for any assumed reionization history, by 
computing the predicted baryon fraction as a function of mass and 
redshift. The predicted M c (z) evolution is in very good agreement 
with the full simulations. 

Our simulations and models both neglect self-shielding 
dSusa&Umemu^ l2004ah. and radiative transfer and non- 
equilibrium effects Abel & Haehneltl 1999). We have performed 
a numerical experiment in which we artificially increase the 
photoheating rate at higher density, to estimate how our results 
might be changed by these effects. The modified heating rate does 
not change the mean IGM temperature and hence also the filtering 
mass remains the same, yet M c determined from the simulations 
increases significantly. Our models correctly predict this increase, 
demonstrating that the key physical quantity regulating the baryon 
fraction in haloes i s neith er the mean temperature of the IGM 
favoured bv lGnedir] (12 000). nor temperature of gas of A ~ 1000 
favoured by Hoeft et all but the temperature of the gas hanging 
around haloes (A ~ l/3A v j r ). 

Our models are easy to implement in any semianalytic model 
of galaxy formation, provide one can compute the equilibrium tem- 
perature at A = l/3A v ir as a function of redshift for the assumed 
UV-background. 
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APPENDIX A: AN IMPROVED MODEL FOR GAS 
ACCRETION ONTO HALOES 

Once density of accreting gas, p acc , is specified, its temperature can 
be computed as T acc = r oq (p acc ). In model A, we assumed that 
the density of gas accreting onto a halo is one-third of the virial gas 
density, i.e. p acc = l/3p v ir(z). If the temperature at this density , 
T e q(pvir/3), is higher than the halo's virial temperature, T v i r , the 
gas cannot accrete. If merging haloes have an unusually low baryon 
fraction, it may be that a larger fraction of baryons is available for 
accretion. We attempt to include this idea in our model 'B' , by 
writing the density of accreting gas as 



- a Macc 
p a cc — p MmaK Pvir, 



(Al) 



where (3 is a free parameter of the order of unity. A halo with baryon 
poor progenitors has higher M acc /M™ ax , and hence higher p acc . 
The temperature of the accreting gas is now evaluated at p acc , and 
can be higher than the value p v i r /3 of model A. A value of (3 — 
2/3 improves the predictions of model B over that of the simpler 
model A. With this value of j3, models A and B accrete gas at the 
same density for a halo with mass M c , which has a baryon fraction 
(/b)/2. 



APPENDIX B: FILTERING MASS 

iGnedin & Huil |l998) found that growth of density fluctuations in 
the gas is suppressed for comoving wavenumber k > /cf, where 
the critical wavenumber fcp is related to the Jeans wavenumber fcj 
by 



,D(t')+2H(t')D(t') 



dt" 



(Bl) 



fcf(i) D(t) J k*(t>) U « 2 (i") 

The physical Jeans wavenumber fcj is defined as 

kj = -(4^G(p tot }) 1/2 . (B2) 

Here, D(t) and H(t) are the linear growth factor, and Hubble con- 
stant as functions of cosmic time t, respectively, and the sound 
speed c s is defined as 

t^T. (B3) 
^3 (im p ) 

Note that lGnedinl d2000h used the volume-averaged IGM tempera- 
ture, {T)v, instead of To. He compared the corresponding filtering 
mass 

4tt . , /27rcA 3 



m f = - {ptot ) I : 



(B4) 



10 T. Okamoto, L. Gao and T. Theuns 




Figure Bl. Upper panel: Evolution of the IGM temperature at mean den- 
sity (solid line) and the volume-averaged temperatures in the 'reference' 
(dotted line) and 'low-resolution' (dashed line) simulations. These temper- 
atures are used to compute the filtering mass shown in the lower panel. 
Lower panel : Comparison of filtering mass, Afp , obtained from Eq. )B4K 
and characteristic mass, M c , obtained from the simulations. The filtering 
mass Mp « 10Af c by z = 0, illustrating the fact that the filtering mass 
strongly overestimates the effect of photoheating on the formation of small 
galaxies. 



with the characteristic mass measured in cosmological simulations 
of reionization. He found a good agreement between the character- 
istic mass and the filtering mass down to z = 4. 

In the upper panel of Fig. IB II we show evolution of To and 
the volume-averaged temperatures in the 'reference' and 'low- 
resolution' simulations. Since larger objects form in the 'low- 
resolution' simulation owing to its larger box size, the IGM tem- 
perature in the 'low-resolution' simulation at low redshift becomes 
slightly higher than that in the 'reference' simulation because of the 
shock heating. In the lower panel, we compare the filtering mass 
calculated from both To and (T)v with the characteristic mass. 
Clearly, the filtering mass significantly overestimates the character- 
istic mass at low redshift, while there is a reasonably good agree- 
ment at high redshift (z > 6) between them. 

This paper has been typeset from a IfX/ ET^X file prepared by the 
author. 



